function Hv = H(a,b)
%H возвращает вектор фокусировки для нашей антенной системы 
%   входные аргументы - азимут и угол места


[kx,ky,kz] = sph2cart(b, a, 1);
k = [kx; ky; kz];


r{1} = [1/4, -1/4, 0];
r{2} = [-1/4, -1/4, 0];
r{3} = [-1/4, 1/4, 0];
r{4} = [1/4, 1/4, 0];
NAnt = length(r);

phi= nan(NAnt, 1);

for i= 1:NAnt
    phi(i) = 2*pi*r{i}*k;
end
Hv = exp(1i*phi);
end

